Optical/infrared flares of GRB 080129 from late internal shocks 
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ABSTRACT 



Strong optical and near-infrared (NIR) flares were discovered in the afterglow 
of GRB 080129. Their temporal behaviors, the sudden emergence and the quick 
disappearance, are rather similar to that of many X-ray flares (for instance, the 
^ ! giant flare of GRB 050502B). We argue that the optical/NIR flares following 

GRB 080129 are a low energy analogy of the X-ray flares and the most likely 
interpretation is the "late internal shock model". In this model, both the very 
sharp decline and the very small ratio between the duration and the occurrence 
time of the optical/NIR flares in GRB 080129 can be naturally interpreted. The 
initial Lorentz factor of the flare outflow is found to be ~ 30, consistent with the 
constraint < 120 set by the forward shock afterglow modeling. Other possibilities, 
like the reverse shock emission or the radiation from the continued but weaker 
^ ! and weaker collision between the initial GRB outflow material, are disfavored. 
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1. INTRODUCTION 



GRB 080 129 was triggered a nd located by the Swift Burst Alert Telescope (BAT) at 
06:06:45 UT ( Immler et al. 2008h . The dur ation of prompt emission T 90 is 48 ± 10 sec in 



15 — 350 keV band (IBarthelmy et al. 20081 ). The time-averaged spectrum is best fitted 



by a simple power-law model, whose power law index is 1.34 ± .26. The fluence in the 



15-150 keV band is 8.9±1.4x 10~ 7 erg/cm 2 flBarthelmv et al. 20081 ). The BAT observations 



lasted until 320 seconds after the trigger, then slewed to another location of the sky. X- 
ray telescope (XRT) and the UV-optical telescope(UVOT) started to point to GRB 080129 
until 3.2 x 10 3 s after the trigger. A fading X-ray source was discovered and no emission 
was seen with UVOT. No flare was observed in X-ray band since XRT started to observe 



flHolland 20081 ). 
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The optical/NIR ob servations imaged by GROND started immediately after the trig- 



ger (IGreiner et al. 20081 ). The first images immediately revealed a strongly flaring source. 
Distinguished optical/NIR flares were observed with amplitude ~ 3 mag, duration of 80 s 
(full-width at half maximum; FWHM, hereafter we define the FWHM as the observed vari- 
ability timescale 5t in the flare), peaking at t p ~ 540 s after the GRB trigger. Their rise 
and the decline can be well approximated by t 12 and t~ 8 , respectively. Thereafter, the after- 
glow brightness is continuously rising until 6000 s after the GRB. The optical spectroscopy 
suggests a redshift z = 4.349 for GRB 080129. 

Greiner et al.(2008) interpreted the optical/NIR flares as the radiation of continued 
but weaker and weaker collisions between the material ejected during the prompt emission 
phase. In this work we do not follow their treatment for the following arguments: (i) In such 
a scenario, the NIR/optical flares emerge when the synchrotron self-absorption frequency 
drops below the observer's frequencies. If correct, the NIR and optical flares should have an 
observable/significant time delay, that is the higher the observer's frequency, the earlier the 
arrival time. However, we did not see such a delay in the data (Greiner et al. 2008; see also 
our FiglTJ). (ii) The NIR/optical flares appeared and then peaked at a time t ~ 540s 3> T 90 . 
If the IR/optical flares are indeed from the outflow material ejected during the prompt 
gamma-ray emission phase, their declines are governed by the high latitude emission and 
can not be steeper than (t — Tgo) _( ' 2+/3 ^ ~ t~( 2+l3 \ again inconsistent with the data, where 
(3 < p/2 is the spectral ind ex, y is the power-law index of the energy d istribution of the 



shock-accelerated electrons (IKumar &: Panaitescu 200CH ; iFan fc Wei 20051 ). This puzzle can 



be solved if the jet is so narrow that we have seen its edge, i.e., 9 } < 0.01 (100/ri), where 
9j is the half opening angle and T; is the initial Lorentz factor of the outflow. However such 
a possibility has been convincingly ruled out by the late ti me afterglow observa tion because 



the jet break at ~ 1.8 x 10 4 s suggests a 0j ~ 0.076 > 1/Ti (IGreiner et al. 20081 ). The latter 
argument applies to the reverse shock emission model as well. That's why we won't discuss 
such a possibility in this work, either. 

We note that the temporal behavior of the NIR/optical flares detected in GRB 080129 
is quite similar to that of X-ray flares observed in a good fraction of Swift GRB afterglows 
(e.g., Guetta et al. 2006; Chincarini et al. 200 7). For comparison p urpose, we re-plot both 



the giant X-ray flare following GRB 050502B ([Burrows et al. 20051 ) and the NIR flares in 
GRB 080129 in FigJTJ The physical parameters are summarized in TabJTJand the similarities 
are evident. Motivated by these similarities we suggest that the NIR flares detected in GRB 
080129 should have the same origin of the flares observed in the X-ray a fterglows, i.e., the 



NIR flares should b e powered by the s o-called late internal shocks, too (IFan fc Wei 2005 



Burrows et al. 20051 ; IZhang et al. 20061 ). Such a model was thought to have been ruled 



out by the request of a very large initial Lorentz factor (I^ ~ 800) of the flare outflow 
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( IGreiner et al. 20081). We'll show in this work that R ~ tens, a typical value taken in the X- 
ray flare modeling (jFan &: Wei 20051 ). is large enough to reproduce the data and is consistent 
with the upper limit (< 120) set by the forward shock optical afterglow modelling of GRB 
080129. 



Table 1: The main parameters governing the opical/NIR flares in GRB 080129 and the X-ray 
flare in GRB 050502B. 



flare(s) in GRB 


rise index (an) 


decline index (—012) 


t P (s) 


St/tp 


GRB 080129 


12 


-8 


540 


0.15 


GRB 050502B 


9.5 


-9.0 


740 


0.14 



This work is arranged as the following. In section 2 we briefly introduce the late internal 
shock model and then discuss the identity of two independent constraints that are widely 
used to rule out other possibilities. In section 3 we apply the late internal shock model to 
the IR flares following GRB 080129. We summarize our results with some discussions in 
section 4. 



2. THE LATE INTERNAL SHOCK MODEL 



In the standard fireball model, the GRB prompt emission is powered by the interaction 
of shells with different Lorentz factor s in the relativistic outflow launched by the central 
engine, i.e., the internal shock model (IPaczynski fc Xu 19941 ; iRees fc Meszaros 1994T) while 
the afterglow is believed to be the external forward shock emission (IPiran 19991 ). However, 
since the launch of Swift satellite, energetic X-ray flares have been detected in about half of 
the GRBs afterglows. The temporal behavior of most X-ray flares share some similarities with 
the prompt soft 7— ray emission and can not be interpreted by the external forward shock 
model (See Meszaros 2006; Zhang 2007 for recent reviews). The most likely interpretation is 
the so-called "late internal shocks model" , in which the GRB central engine restarts after the 
prompt emission phase and launches unsteady outflow. The underlying physical processes 
are less clear. Among the various models put forward (see Zhang 2007 for a review) fallback 
accretion onto the nascent black hole may be the most natural one. The collision between 
the fast and slow material of the new outflow can power strong flares peaking in X-ray 
or far-ultraviolet band. The duration of these flares (St) is determined by the re-activity 
process of the central engine and can be much shorter than the occurrence time of the flares. 
On the other hand, since the ejection time (~ £ e j e ) of the last main pulse of the flare is 
close to £ p , the net flux of the high latitude emission of the pulses can be approximated 
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Fig. 1. — The light curves of NIR flares in GRB 080129 compare with that in 
GRB 050502B in X-ray band. Left: The red circles, green pentagons and blue 
points represent observed data from three n ear-infrared channels K s xl0, Hx5, Jx2.5 
on GROND, respectively jGreiner et al. 20081 ). Right: Data get from UNLV GRB 
Group (http://grb. physics. unlv.edu/). The main parameters describing these flares 
are summarized in TabJH 
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by (t p — t C j e ) _ ( 2+/3 \ which can be much steeper than t~( 2 +^) . So the late internal shock 
model can naturally account for the main characters, the sudden emergence and then a 
rapid drop, of the X-ray flares detected so far. For the emissi on of the external sh ocks, 



it is well known that (1) St/t has to be in order of 1 or larger (INakar &: Piran 20031 ); (2) 
the decline can not be steeper than t~( 2+ ^) unless the edge of the GRB ejecta is visible. 
This is because the GRB outflow is curving and emission from high latitude (relative to 
the observer) will rea c h us at later times and giv e rise to a decline shallower than t~( 2+/3 ) 



( iFenimore et al. 19961 : iKumar fc Panaitescu 2000l ). Usually these two limitations have been 



taken as independent evidences for the late internal shock model (e.g., Chincarini et al. 
2007). Below we show that they are highly relevant and even identical)^ 

As shown in Figf2[ f p is the maximum flux at the peaking time t p in the flare, t\ 
and ti are the time at which the flux is half of / p in the rising and decaying light curves, 
respectively. The FWHM time is therefore 8t = t<i — t\. Before and after t p , the light 
curves are approximated by t ai and t~° 2 , respectively. It is straightforward to see that 

= ^faV-^gftif ' and a2 = iSl^FS^ - After some sim P le al g ebraic we have 

6t/t p = 2 1/Q2 - 2- (1/ai) (1) 

Obviously, 5t/t p is irrelevant to f p . For fixing a±, the relationship between St/t p and «2 
can be found in Figj3j with which we can see that the steeper the decay, the smaller the 
St/t p . Particularly for a\ 3> 1, we have a 2 ^> 2 + (3 (hereafter "the decline constraint") and 
St/t p <C 1 simultaneously, suggesting that the two constraints widely used in supporting late 
internal shock model are highly rel evant. This naturally accounts for the fact that many 



X-ray flares satisfy both limitations (jChincarini et al. 20071 ). Please note that our conclusion 



is independent of the underlying physical processes. 

For the purpose of identifying the afterglow emission powered by the central engine, the 
decline constraint may be more g eneral. For examp le, the very sharp drop detected in the 



X-ray afterg low of GRBs 070110 flTroia et al. 20071 ). 060413, 060522, 060607A and 080330 



(IZhang 20091 ) also favors a central engine origin though the constraint St/t <C 1 is violated. 



For the optical/NIR flares of GRB 080129, the peak time is ~ 540 s after the trigger, 
the duration (i.e., the FWHM) is 80 s. We have St/t p ~ 0.15 <C 1 and a 2 = 8, indicating a 
central engine origin of these flare photons. Below we take the late internal shock model to 
reproduce the data. 



1 Fan et al. (2008a) pointed out this in a proceeding paper but did not prove it. 
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Fig. 2. — A schematic plot of a flare in the GRB afterglow. 
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Fig. 3. — The relationship of decay power law index a<i and the ratio of 5t/t p . The rising 
power law index ot\ is taken as 2, 4, 6, 8, 10, 12, 14, respectively. The horizontal dot line 
represents = 2 + (3 = 3.25, where we take /3=p/2=1.25, the vertical dot line represents 
5t/t p = 1. 
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3. Physical parameters and the synchrotron radiation of GRB 080129 in late internal 

shock model 

We assume that the Lorentz factors of the ejected material in the re-starting outflow are 
highly variable, and take T s ~ 10 and Tf ~ 100 as the typical Lorentz factor of the slow and 
fast shells, respectively. The masses of the fast and slow shells are taken as nif ~ m s . In the 
late internal shock model, the inner fast shell will catch up with the outer slow shell at the 
radius ~ 2r 2 c<5ti/(l + z) (where <5ti is taken as the observed typical variability timescale of 
one pulse in GRB 080129 optical/NIR fla re), and inte rnal shock are generated. The merged 



shell's Lorentz factor is Ti « V / TfTs ~ 30 (jPiran 19991 ). and the Lorentz factor of the internal 



shock can be estimated as 7 s h ~ (yTfTT^ + ^/r s /Tf)/2. 

Adopting the cosmological parameters H = 70kms~ 1 Mpc~ 1 , Q M = 0.3, and Qa = 0.7, 
we have a luminosity distance Dl = 1.2 x 10 29 cm for GRB 080129 at a redshift z = 4.35. The 
observed maximum flu x of the flare in GRB 080129 is about~ 0.5 mJy (NIR) and ~0.3 mJy 



(optical), respectively (IGreiner et al. 20081 ). Assuming an efficiency factor of the optical flare 
e ~ 0.1, the total luminosity of the flare outflow ican be estimated by L m ~ 3 x 10 48 ergs -1 . 
This luminosity implies that the fallback accretion rate is about ~ 10~ 5 — 10~ 3 times that 
of the GRB prompt accretion, if the efficiency factor of converting the accretion energy int o 



the kinetic energy of the outflow is nearly a constant (jMacFadyen. Woosley fc Herger 200ll ). 



The variability timescale 8t\ of GRB080129's optical flare is significantly longer than 
that of the prompt emission. Here we take 5t\ ~ 30s, as suggested by the smoothness of the 
flare light curves. The typical radius of the late internal shock is R int ~ 2Tfc5ti/(l + z) ~ 
3.4 x 10 14 Tf l 5 8t- h i_5 cm. In this work we take the convenience Q x = Q/10 x in units of cgs 
unless with specific notation. Below, following Fan & Wei (2005), we show that with the 
parameters F{ ~ 30, L m ~ 3 x 10 48 erg s _1 , e e = 0.4, e-Q = 0.05, and 8t\ = 30s, the flare data 
can be reasonably reproduced. 

We first investigate the interaction between the inner fast shell and the outer slower 
one. The comoving number density of the electrons is n e ~ L m / (4nrf R 2 m p c 3 ) ~ 4.6 x 
10 7 L mi 4 8 . 5 rf 1 6 5 5t 1 ~ 1 2 5 , where m p is the rest ma ss of proton. The thermal e nergy density of the 



shocked material is e = 47 s h(7 s h — l)n e m p c 2 (IBlandford fc McKee 19761 ). So the strength of 
magnetic field can be estimated as 

B » (87r eB e) 1 / 2 ^6.7xl0 2 G e ^ L3 (^) 1 / 2 

( 7sh - l^Zs.^J^L- (2) 

As usual, we assume that in the shock front, the accelerated electrons take an energy 
distribution dn e /dj e oc % p for j e > y e>m , where 7 e>m = e e ( 7sh - l)[(p - 2)m p /(p - l)m c ] is 
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the minimum Lorentz factor of the shocked electrons (jSari et al. 19981 ). and m e is the rest 
mass of an electron. Here for GRB 080129, we take p = 2.5. We can get the observed typical 
frequency of the synchrotron radiation 



v m = jl m q e TB/[2(l + z)-nm e c\ 

3.7 x 10 14 Hz e e 2 _ . 4 6^_ L 3( 7sh - l) 5 / 2 ( 7sh /2) 1 / 2 

^m,4S.5^i,1.5^h,l.57 (3) 

where q e is the charge of the electron. 

The cooling Lor entz factor is es timated by 7e, c ~ 7.7 x 10 8 (1 + z)/(TB 2 5ti). So the 



cooling frequency is (jSari et al. 19981 ) 



7 2 g c rB/[2(l + z)vrm e c] 



~ 1.0 x 10 12 Hz e^LC— )" 3/2 

(7sh - ir^L-JilTlJt-^. (4) 



The synchrotron self-absorption frequency can be estimated as 



u a » 1.2 x 10 14 Hz e^L^lsKTsh 

-8/7,-5/7 
1 i,1.5 U H,1.5 ■ 



i \ T^hil/14 

1 2 J 



(5) 



The maximum spectral flux of the synchrotron radiation is -F max m 3\^3& p (l+z)N e m e c 2 aTTB /(32ir 2 q e L 
where N e is the total number of emitting electrons, N e = L m 8t\/[(1 + z)T\m v c 2 } = 9.3 x 
10 50 L nl) 48.5rr' 1 _ s 5ti j i.5, where 5U is the observed typical variabil ity timescale of the tota l flare, 
and $ p is a function of p. For p = 2. 5 we have Q p = 0.6 (IWijers &: Galama 19991 ). For 
v c < v a < v < u m , the predicted flux is (jSari et al. 19981 ) 



-1/2 



5.0 x 10~ 4 Jy [i//(3.0 x 10 14 Hz )]~ 1/2 e 

(7 sh - i)- 1/4 L 3 / 4 4 8 . 5 5t?/ 1 2 5 r iiL5 D 



-1/2^1/4 ^Tsh^-1/4 



B.-1.3' 



-2 
L,29 



(6) 



Taking z^nir = 3.0 x 10 Hz , we have F Umn 
GRB 080129's flare in near-infrared band. 



0.5 mJy, consistent with the observation of 
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In the optical and X-ray band satisfying v c < z/ a < v m < the flux can be estimated as 
F„ = F max (z/ m /z/ c )- 1 / 2 (z//z/ m )~P/ 2 JSari et al. 1998h . Taking z/ opt = 5.0 x 10 14 Hz for optical 
band and 2 x 10 17 Hz for X-ray band, we have F u . ~ 0.3mJy and F Ux ~ 1.7 x 10 _4 mJy, 
respectively. Approximately the optical peak flux of the flare is ~ 0.3mJy, as inferred from 
Fig.l of Greiner et al.(2008). So our result is consistent with the optical data, too. In the 
X-ray band, no observation was carried out for t < 3.2 x 10 3 s. So it is impossible to test 
our predication in X-ray band. 



4. DISCUSSION 

In ~ 10 2 — 10 5 s after the trigger of GRBs, bright X-ray flares have been well detected 
in a good fraction of Swift GRB X-ray afterglows (Falcone et al. 2007; Chincarini et al. 
2007). However, for many X-ray flares the peak energy is unknown and the upper limit is 
about 0.2 keV. Fan & Piran (2006) speculated that some X-ray flares actually peaked in 
UV/optical band and thus should be classified as UV/optical flares. However, before 2008 
people had not detected a canonical optical flare with plenty of data. The best candidate 
of UV/optical flare may be that detected in GRB 050904 (Boer et al. 2006), for which, 
unfortunately, the reverse shock model can not be ruled out (Wei, Yan & Fan 2006). The 
situation cha nged dramatically af ter the release of the early optical/NIR afterglow data of 



GRB 080129 flGreiner et al. 20081 ) 



The optical/NIR flares following GRB 080129 have very sharp decline («2 3> 2 + (3) 
and very small 5t/t p (~ 0.15), rather similar to that of the giant X-ray flare following GRB 
050502B. These two characters rule out the possibility of being the reverse shock emission 
or being the radiation of the continued but weaker and weaker collision between the outflow 
material ejected during the prompt emission phase. Instead, these optical/NIR flares can 
be attributed to the re-activity of the central engine, as the X-ray flares detected in a good 
fraction of Swift GRB X-ray afterglows. In the framework of late internal shock model, 
with reasonable physical parameters (in particular Ti ~ tens) we calculate the synchrotron 
radiation. The typical frequency is just in near infrared band and the flux estimated in 
near-infrared and optical band are also consistent with the observations (see section 3 for 
details). We conclude that the flares in GRB 080129 peaking in NIR/optical band are a low 
energy analogy of the X-ray flares, confirming the speculation of Fan & Piran (2006). 

The identification of a low energy analogy of X-ray flares in optical/IR band also helps 
the people to diagnose the physical composition of the outflow launched by the re-activity of 
the central engine. Fan et al. (2008b) showed that polarimetry of the flares is highly needed 
to achieve such a goal. Technically the optical polarimetry is much more plausible than the 
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X-ray polarimetry at present (jCovino et al. 19991 ). 

In this work we also show that the two constraints a 2 3> 2+ (3 (i.e., the decline constraint) 
and St/t <C 1, widely /separately used to support the "central engine origin" of the afterglow 
emission, are highly relevant and even identical (see section 2 for details) for the flares. The 
decline constraint may be more g eneral. For examp le, the very sharp drop detected in the 
X-ray aftergl ow of GRBs 070110 faroia et al. 2007h . 060413, 060522, 060607A and 080330 
( iZhang 20091 ) is in support of a central engine origin though the constraint 5t/t <C 1 is 
unsatisfied. 
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